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ABSTRACT 

The B2Vn star HR7355 is found to be a He-rich magnetic star. Spectropolarimetric data were 
obtained with FORS1 at UT2 on Paranal observatory to measure the disk-averaged longi- 
tudinal magnetic field at various phases of the presumed 0.52 d cycle. A variable magnetic 
field with strengths between (B z ) = —2200 and +3200G was found, with confidence limits 
of 100 to 130G. The field topology is that of an oblique dipole, while the star itself is seen 
about equator-on. In the intensity spectra the Hel-lines show the typical equivalent width vari- 
ability of He-strong stars, usually attributed to surface abundance spots. The amplitudes of 
the equivalent width variability of the Hei lines are extraordinarily strong compared to other 
cases. These results not only put HR7355 unambiguously among the early-type magnetic 
stars, but confirm its outstanding nature: With wsini = 320 km s -1 the parameter space in 
which He-strong stars are known to exist has doubled in terms of rotational velocity. 
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1 INTRODUCTION 

iRivinius et all < l2008h suggested the B2Vn star HR7355 
(HD 182 180, HIP 95 408) to be a He-rich magnetic star 
with a magnetosphere containing trapped gas tha t pro- 
duces hydrogen line emission ([Townsend & Owockil I2005L 
ITownsend. Owocki & Ud-Doulal l2007h . They based their sugges- 
tion on evidence derived from two FEROS echelle spectra and a 
Hipparcos light-curve, with a period of either single-wave 0.26 d 
or double-wave 0.52 d. However, a period of 0.26 could not be 
due to rotational modulation, since this would require a rotational 
speed well above the critical threshold. 

If confirmed, HR7355 would have a unique position in the 
He-st rong class: With v sin i = 320 km s" 1 dAbt Levato & Grossd 
l2002h it would be the the most rapidly rotating magnetic star in the 
upper HR-diagram; by a factor of about two ahead of t he current 
record-holder, a Ori E, which has v sin i = 165 kms" 1 dAbt et all 
op. cit.). This would make the star a show-case for the Rigidly Ro- 
tating Magnetosphere model, that so far, t hough with great suc- 
cess, was applied to only one star, a Ori E ( ITownsend & Owocki 
op. cit.). In order to test this claim, a spectropolarimetric campaign 
was carried out in 2008. 



* Based on observations made with ESO Telescopes at the Paranal Obser- 
vatory under programme ID: 081.D-2005 
f triviniu@eso.org 



2 OBSERVATIONS AND DATA REDUCTION 



The data described below were obtained a t the VLT with the 
FORS1 instrument t Appenzeller et al.lll998h . which is equipped 
with a super-achromatic quarter-wave retarder plate. The instru- 
ment setup used the 1200B holographic prism with a dispersion 
of 24 A/mm over a range of 366 to 51 1 nm. At a slit-width of 0.3", 
a (binned) pixel scale of 0.25"per pixel, and a (binned) pixel size of 
30 fim, this gives about 0.9 A per pixel, which slightly oversamples 
the resolution element of about 1.3 A. 

Each individual measurement (see TablefJJ consisted of a se- 
quence of eight exposures, taken at retarder plate position angles in 
the sequence (-45° , -45° , 45° , 45° , -45° , -45° , 45° , 45°) with 
respect to the axis of the Wollaston prism. The exposure time was 
two seconds for the respective 8 exposures. The data was taken 
during four nights in service mode. During two of these nights 
the sequence was repeated, so that we have obtained six data 
sets, each one reaching an approximate S /TV-ratio of about 1000 
per resolution element. The method for extracting the Stokes V 

[ arameter from the da ta is de scribed in the FOR S user manual 
lehin & O'Brienll2008l) and bv lBagnulo eTal] d2002l) . 

Example Stokes V data are presented in Fig.fJJ in this figure, 
and for the field strength determination (as in Fig. |2), there is no 
rebinning from pixel space to uniform wavelength space. We cal- 
culated the wavelength for the respective data points based on the 
HgCd and He arc line spectra taken during daytime. 

Flat fielding is largely unnecessary for the Stokes V determi- 
nations, because the alternate measurement of —V/I and V/I at 
—45° and 45° retarder plate angles cancels out any pixel-to-pixel 
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Figure 1. Normalized circular polarization V/I of HR7355 at 
MJD=54669.33, for which a longitudinal magnetic field component of 
3200 G was derived. Top: The spectra around the ten Balmer lines from 
H/3 to H13. Bottom: The V/I vs. the bracketed part of the right term in 
Eqn.[T]and the linear regression to the data as a solid line. The typical error 
per measurement is shown as a bar; the error in the abscissa is less than 
the symbol size. The mass of points clustering around (0,0) are due to the 
unpolarized continuum. 



variations in the CCD response; nevertheless, because intensity (I) 
data are required to derive the field strength, we applied an appro- 
priate flat-field correction to all spectra. 



3 RESULTS 

3.1 Magnetic field 

As discussed by Bagnulo et al. ( 2002), under the weak field ap- 
proximation one can derive the mean longitudinal magnetic field 
component over the stellar disk (B z ) from the circular polarimetry 
and the gradient of the intensity spectrum with the equation: 

2 dl 



-KeffCzA^ 



JdA 



x (B z 



(1) 



Here we used a factor C z summarizing the physical constants 
like the electron charge e, the electron mass m c and the speed of 
the light c in the form: 



C, 



47rm c c 2 



4.6710" 



and an effective Lande factor g c ff = 
ICasini & Landi D egl'Inno centil jl994f) . 



(2) 

1 as discussed by 




Figure 2. Minimum of the two-dimensional probability density disctribu- 
tion computed from the artificial test data (see Sect. |3.U for the observed 
data shown in Fig.[T] and an overplot of the x 2 -contrours enclosing 68.3%, 
95.4% and 99.7% probabilities (upper panel) and one dimensional proba- 
bility distribution of the derived field strength (lower panel). 



It should be clearly mentioned that there are a number of ap- 
proximations and assumptions in equation [T] above. One of those 
is the week field approximation for the magnetism, i.e. a Zeeman 
split of less then the intrinsic line thermal and pressure broadening, 
which is, however, safely applicable to a field of only a few kilo- 
gauss in a main sequence B star. A concern prompting far more cau- 
tion when interpreting the derived numbers is that the spectra are, in 
fact, an average over the stellar disk longitudinal fie ld component, 
witho ut considering e.g. the limb darkening (See Bag nulo et al.l 
for a full discussion of the limitations and caveats of the 
method). However, as this quantity, labelled (B z ), is typially the 
most easily derived m agnetic obse rvable, it is also the most pub- 
lished in the literature I Wade 20031) . and as such easily comparable 
to other measurements. A more complete magnetic modelling to 
derive the physical dipole field strength is left to a later work. 

We used a % 2 minimization approach to fit the linear model 
V/I — A + (B z ) x x to the observational data; this model is 
based on eqn. Q] with the addition of a constant term A. The ad- 
ditional term A is physically expected to be zero, as it is the V / 1 
of the continuum. However, the employed measurement principle 
does not necessarily guarantee A = in the measured data, so we 
allow non-zero values in order to improve the linear regression to 
the slope. It actually turns out that the derived values for A are con- 
sistently above zero by between 0.4 x 10 -4 and 2.5 x 10~ 4 with 
a typical a — 0.5 x 10~ 4 , and on average 1.5 x 10~ 4 . This likely 
is an instrumental effect. In any case this is rather small compared 
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MJD (B z ) a Bz W x , Hel 4388 Wx, Hel 4713 phase 

[G] [G] [mA] [mA] 



54652.327 


-2210 


± 130 


1720 


530 


0.681 


54656.078 


-1760 


± 130 


850 


310 


0.875 


54656.146 


230 


± 100 


740 


280 


0.006 


54661.327 


-1110 


± 120 


660 


250 


0.942 


54669.186 


220 


± 110 


820 


310 


0.013 


54669.330 


3200 


± 130 


1530 


480 


0.289 



Table 1. List of magnetic field and equivalent width measurements derived from the FORS1 data. The typical error for the equivalent width is about 5 %. For 
computing the phase we use the ephemeris defined in Eqn.[5] Only Balmer lines were included in the field measurements. 



to the V/ 1 peak-to-peak amplitude induced by the magnetic field, 
which is up to 40 times higher (see Fig.[TJ.. 

The formal errors of V/I, which are derived from the well 
known photon- and detector noise characteristics, were used to as- 
sign weights for the \ 2 fit- The derived (B z ) are given in TableQ] 

Given recent concerns expressed over the error b ars associated 
with FORS1 magnetic field measurements (see, e.g. lSilvester et al.l 
2009), we have exercised particular care in deriving the errors listed 
in Table[T] Applying the bootstrap Monte Carlo approach described 
bv lPressetatl fl992), we used the observational data from each 
epoch to construct one million corresponding synthetic datasets. 
We then applied the same fitting approach as above to obtain A and 
(B z ) values for every synthetic dataset. These values exhibit a dis- 
tribution about those derived from the actual observations; Fig. [2] 
shows this distribution for the MJD 54 669.330 observations, both 
as a 2-dimensional probability density map, and as a 1 -dimensional 
probability density function in (B z ). Plotted over the map are con- 
tours of constant \ 2 that enclose 68.3%, 95.4% and 99.7% of the 
synthetic datasets; likewise, plotted over the probability density 
function are the symmetric bounds enclosing 68.3%, 95.4% and 
99.7% of the synthetic datasets. The bounds for the 68.3% case are 
assigned as the error bars quoted in Table [T] because the proba- 
bility density function in Fig.|2]is close to Gaussian (and likewise 
for the data from other epochs), these error bars can be regarded 
as approximate 1-cr confidence intervals. Note, however, that we 
have not made any specific assumptions regarding the propagation 
of errors through our modeling process; the Monte Carlo approach 
naturally results in error estimates that directly reflect the combined 
characteristics of the observations and the modeling. 



3.2 Spectral lines 

Next to the Stokes V spectra, we examined the intensity spectra 
observed by FORS1. Due to the very high S/N and the large rota- 
tional velocity, even the relatively low dispersion of FORS 1 allows 
high quality equivalent width (EW) measurements of the spectral 
lines. 

The Hi lines, with the exception of H/3, show little variabil- 
ity, except in the very core. The spectra taken at MJD=54656.146 
and 54669. 1 86 show enhanced core absorption. At these epochs the 
(B z ) had null-detections, which is in agreement with the crossing 
of the magnetic equator through the line of sight, because this is the 
time at which a corotating cloud of circumstellar matter, located at 
the crossing of the magnetic and rotational equators, is expected to 
pass in front of the star. 

In H/3 this circumstellar absorption is seen as well, but in ad- 
dition the other four spectra show variability in the line wings, in 
particular in the red one. This is probably a signature of the cir- 
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Figure 3. Intensity spectra measured by FORS. The Hel-variability is well 
seen, while the Balmer lines vary only little and the CII4267 and Sim 4552 
lines remain stable. 



cumstellar emission ari sing in the corotatin g clouds, well seen in 
the two FEROS spectra dRivinius et al.ll2008h for Ha. 

All Hel lines in the observed range follow a similar variation 
pattern, varying both in strength and in profile. The lines show- 
ing strong broadening wings have a larger amplitude in Wx due to 
variation in these wings, but the profile variability is better seen in 
weaker lines, like Hel 47 13. 



3.3 Forbidden Hel lines 

In addition to the well known stellar Hi and Hel lines, there is 
significant variability at A ~ 4045 A and bluewards of Hel4922, 
at A ~ 4911 A. We identify the features with the forbidden lines 
Hei4045 (2 3 P 5 3 P) and Hei 4911 (2 1 P -> 4 X P), features w ell 
known in extreme Helium stars (Beau champ & Wesemaelll 998). 

In the FORS data, variability due to forbidden components of 
Hel at AA 4045 and 4911 was found. At a closer inspection, this 
spectral signature is also present in the old FEROS spectra, and a 
search in other magnetic He-strong stars reveals the presence of 
these lines also in a Ori E and V 1046 Ori. Less certain, though 
not excluded from our archival data, are those lines in HD 64 740 
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Figure 4. Forbidden Hel 4045,4911 signatures in some known magnetic 
He-strong stars. The forbidden components are clearly seen in in the slowly 
rotating HD 94 446, but also in HR 7355, V 1046 Ori and cj Ori E, and possi- 
bly in HD 64 740 and HD 37 776 as well. Note that these are single spectra, 
so shifts occur as the features are Doppler shifted due to the rotational vari- 
ation. 



and HD 37 776 (Fig. [5}. These lines are never seen in B stars with 
normal He abundance. 

In other words, the presence of these lines indicates a Helium 
overabundance wrt. solar values. Such He-strong stars, especially 
when having rotationally broadened lines, are quite hard to diag- 
nose without detailed abundance analysis, as has also happened 
for HR7355, instead of classifying it as chemically peculiar it was 
rather classified as a later type that it actually has. This means the 
presence of the two forbidden Helium lines can safely be taken as 
an indicator for a He-strong star. As many He-strong stars are mag- 
netic (and He-weak stars are typically of later spectral type), with 
the same reasoning we consider the presence of these lines as an 
indicator for a magnetic field in early B-type stars. 



4 DISCUSSION 
4.1 Ephemeris 

As we know from a Ori E teeiners et alj200(jh . the phase most easy 
to pin down, both spectroscopically and photometrically, is the cen- 
ter of one of the two circumstellar clouds crossing in front of the 
star. We chose this point also for HR7355 as (f> — 0. A correspond- 
ing epoch To can then be identified as the time of strongest absorp- 
tion in Ha, i.e. the maximal equivalent width. While our FORS 



data do not cover Ha, and the FEROS data were not taken in this 
phase, a spectroscopic campaign with the echelle instrument UVES 
at Paranal that has been completed in the second half of 2009, sug- 
gests MJD= 54 940.33 as epoch. However, since a full analysis 
of the UVES data is far beyond the scope of this discovery report, 
and the epoch is not as critical as the period for our purpose, we 
leave a detailed analysis and discussion of the Ha equivalent width 
curve to a later work (Rivinius et al, in prep.), and here just adopt 
the epoch. Nevertheless it is reassuring that additional absorption in 
H/3 in the FORS1 data is observed at </> = 0.006 and 0.013, when 
the data are sorted with the period derived below. These two occur- 
rences, 25 cycles apart from each other and more than 500 cycles 
apart from the selected epoch, are in full agreement with this epoch 
derived from UVES data. 

For the period, iRivinius et aH (2008) save a value of 
0.521428(6) d for the Hipparcos data, under the condition that the 
variations were double-wave sinsoidal. The single wave period, 
which in any case would have been to short to be rotational, is 
firmly excluded by the new magnetic measurements. 

However, also the 0.52 d period is not able to satisfactorily 
phase all available data, and so its value had to be improved. In or- 
der to do so, we demand a period to sort all three data types, i.e. 
the photometric data (1990-1993, double wave), the magnetic data 
(2008, single wave), and the equivalent widths (1999-2008, possi- 
bly double wave). With respect to the originally published value, 
the closest period soring all three data-types is 0.521442(4) d. Al- 
though the EW data suffers strong seasonal aliasing, already its 
most nearby alias is excluded by the almost completely scrambled 
photometric curve for this period. The most nearby period value 
for which EW and photometry could be reconciled is incompati- 
ble with the magnetic phase-curve. We thus conclude that the true 
rotational period of HR 7355 must be within the above value's un- 
certainties. 

In a final step, we can assume that the photometric minima do 
have a certain phase relation to the spectroscopic curve. There are 
two possiblities: First, the minima could be due to cloud eclipses 
in front of the star, then we can further require that one of the pho- 
tometric minima occurs at phase cj> — 0. Second, the photometric 
might be due to photospheric flux modulation in the He-enriched 
parts dKrticka et ai]|2007t iMikulasek etall feoiO). In this case, the 
photometric minima would coincide with minimal Hel equivalent 
width. However, in this particular star, the maximal Ha equivalent 
width and the minimal Hel equivalent width are almost simulta- 
neous, so that a final decision about this can only be made with 
new photometry more simultanous with recent observations, in any 
case, already with the period derived from the equivalent widths 
alone one of the two photometric minima is very close to <f> — 0, 
so that under the assumption that the photometric variability is due 
to the circumstellar material obscuring the line of sight the period 
becomes P — 0.521444(3) d. We will use the latter value for the 
discussion, but note that this relies on the assumption that the pho- 
tometric minima are due to eclipses, while the best period without 
this assumption is P = 0.521442(4) d 

The ephemeris used thus is 

T m max. absorption (MJD) = 54940.33 + 0.521444(3) x E (3) 



4.2 Periodic variations 

With the above choice of epoch we expect to see photometric min- 
ima, as well as B z — 0, at phase = 0, i.e. when the mag- 
netic equator is facing towards us (Fig. [5}. Then using the above 
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Figure 5. Measured quantities (see Table[T] also for typical errors) as func- 
tion of the rotational phase (Eqn[5J- The Hipparcos double wave light-curve 
is shown in the uppermost panel, the Hel4388 and 4713 equivalent widths 
below as triangles. The W\ measured in the FEROS spectra IRivin ius et alj 
2008) are shown as red (1999 data) and blue (2004 data) squares. The low- 
ermost panels shows the (B z ) field measurements and the curve defined by 
Eqn.|4] 



a factor of 2 for lines like Hel4713, and even a factor of 3 for strong 
lines with significant broadening wings like Hel4388. This strong 
modulation is indicative for two large Helium enhanced patches on 
the surface close the equator at opposite longitudes, which point to 
a large angle /3 between the rotational and magnetic axes. 



5 SUMMARY 

The results confirm the magnetic and He-strong nature of HR 7355. 
We have detected a magnetic field of multi-kilogauss strength, 
varying with the rotational period, which is P = 0.521444(3) d 
under the assumption of the photometric minima being eclipses, 
and P — 0.521442(4) d without this assumption. In order to avoid 
a potential underestimation of the confidence limits of the magnetic 
fie ld measurement (as s uspected for previous FORS measurements 
bv lSilvester et alj2009h , we derived them with a boot-strap Monte- 
Carlo method, which does not implicitely assume a statistics for 
the error propagation, but numerically reconstructs the probability 
distribution from which the actual observation was drawn. 

The observed magnetic field and the suggested topology 
(see Sect. |4.2b is i n agreement with the initial hypothesis by 
IRivinius et al.1 ( 120081) , namely that of an equatorially seen star, i.e. 
i ~ 90°, with an oblique magnetic dipole. This topology gives 
rise to a double-wave light curve, either as the corotating magne- 
tospheric clouds, magnetically bound at the crossings of magnetic 
and rotational eq uator, sweep thr ough the line of sight twice per 
rotational period |Townsendll2008h . or due to the modulation of th e 
photospheric flux by the abundance pattern ( Mikulasek etalj201Ch . 



ephemeris, and assuming a sinusoidal variation of the magnetic 
field, we estimate the field curve as 

(B x ) (t) = 350 G + 2850 G x sin ^2tt ( - - 0.02^ (4) 

where t is the date, and T and P are epoch and period from Eqn.[3] 
and the shift of 0.02 in phase is required to fulfill the above condi- 
tion B Zi< /, = o = G, due to the constant term of 350 G. 

The f act that -B max fa — B m i n confirms the assumption of 
i ~ 90° bv lRivinius et all J2008h . while the term of 350 G is easily 
explained either by a slight offset from this value, or by an off- 
center magnetic dipole. 

Such an off-center dipole would cause a non-sinusoidal field 
curve. In fact, due to two datapoints being taken at identical phases, 
there are effectively only five points to constrain the three free pa- 
rameters of Eqn. [4] The field curve, although it seems to fit very 
well, is thus not well constrained, which is why we cannot really 
give confidence limits for the parameters unless further data has 
been obtained. In particular, we stress that the exactly sinusoidal 
shape of the curve is rather an assumption than an actually observed 
property. 

When the magnetic poles face the observer, (B z ) becomes 
maximal. Although there are only few points, not sampling the 
curve in all detail and in particualr not necessarily the respective 
maxima and minima, it is clear that the Hel absorption is much 
stronger in these phases than when the magnetic equator is visible 
(Fig. |5j. This is in full agreemen t with the behavior observed in 
other He-strong stars, like a Ori E JReiners e t al. 2000). 

The amplitude of the Hel EW variations is considerably larger 
than in a Ori E, however. There the maximal EW is only about a 
factor of 1.3 to 1.5 stronger than the minimal one, depending on the 
spectral line. In HR7355 the lines strengthen, wrt. their minima, by 
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